Existence of Quasi-stationary states at the 
Long Range threshold 
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Abstract: In this paper the hfetime of quasi-stationary states (QSS) in the a— HMF model 
are investigated at the long range threshold {a = 1). It is found that QSS exist and have a 
diverging lifetime t{N) with system size which scales as T(N)~log7V, which contrast to the 
exhibited power law for a < 1 and the observed finite lifetime for a > 1. Another feature of 
the long range nature of the system beyond the threshold ( a > 1) namely a phase transition is 
displayed for a = 1.5. The definition of a long range system is as well discussed. 
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1. INTRODUCTION 



Studying the dynamics of Hamiltonian systems with a 
• large number of degrees of freedom and its connection to 
' equilibrium statistical mechanics has been a long standing 
problem. The relaxation to statistical equilibrium has been 
' under scrutiny ever since the pioneering work of Fermi and 
I the FPU problem(l). Moreover, since the advent of pow- 
erful computers and for specific systems within a class of 
initial conditions, integrating numerically Hamiltonian dy- 
namics has proven to be competitive in regards to Monte- 
> Carlo schemes for the study of statistical properties (see 
' for instance (2; 3) and references therein). The assumption 
made is that since the system admits only a few conserved 
quantities for generic initial conditions, once the dimen- 
I sions of phase space are large enough, microscopic Hamil- 
, tonian chaos should be at play and be sufficiently strong to 
provide the foundation for the statistical approach within 
the micro-canonical ensemble. However recent studies have 
shown that there is an increase of regularity with system 
size in the microscopic dynamics when considering sys- 
tems with long range interactions (4; 5; 6; 7). Indeed, the 
statistical and dynamical properties of these systems are 
still under debate. For instance extensivity is not always 
provided and discrepancies between canonical and micro- 
canonical ensembles can be found such as negative specific 
heats for the latter (8; 9). Moreover, phase transitions for 
systems embedded in one dimension can be found. 

In particular, long range systems often display a slow re- 
laxation to equilibrium. Starting from an initial condition 
they are in fact trapped in long-lasting out of equilibrium 
regimes, termed in the literature Quasi Stationary States 



(QSS) which have distinct macroscopic characteristics, 
when compared to the equilibrium configuration. 

A now paradigmatic model of long range interactions 
Hamiltonian systems is the Hamiltonian Mean Field 
(HMF) model (10), which corresponds to a mean field XY- 
model with a kinetic energy term (rotators). In the limit of 
infinite system size the HMF model can be described using 
a Vlasov equation (11; 9). More recently, stationary states 
have been constructed using invariant measures of systems 
composed of uncoupled pendula (6), more specifically it 
was emphasized that the microscopic dynamics in the mag- 
netized stationary state is regular and explicitly known. 
This observation lead to explain the abundance of regular 
orbits as revealed in (5). These results first obtained for the 
HMF have been extended for the case when the coupling 
constant depends on the distance between sites, namely 
for a— HMF model in its long range version (a < 1)(7). 
This model was introduced for instance in (12) and and 
displays identical equilibrium features as the HMF(13; 9). 
In fact it was shown that all stationary states of the 
HMF model are as well stationary states of the a— HMF 
model, that microscopic dynamics is as well regular, at the 
price of microscopic spatial complexity, which is locally 
scale invariant (7). Before going on we write the governing 
Hamiltonian of the model: 
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where qi stands for some spin angle located on the lattice 
site z, and pi is its canonically conjugate momentum. The 
distance ||i — j|| is actually the shortest distance on the 



circle of perimeter A'' — 1, so that the systems can be 
isolated and still translational invariant along the lattice. 
The mean field model is recovered for a = 0, and for N 
even, we write 

N/2-1 

iV= ( ^1 +2 V -, (2) 
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to insure extensivity. The equations of motions of element 
i are derived from the Hamiltonian (1): 

Pi = - sm{qi)Ci + cos{qi)St = Ah sin(5i - ipi) , (3) 
ii^Pi, (4) 
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Ci and Si identify the two components of a magnetization 
per site, with modulus Mi 



V^ 



and phase 



tpi = a,Tcta,n{Si/Ci). For large N , and assuming < a < 1, 
we have 

-^(iV/2)i-". (7) 
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We then can use (7) in Eq.(5) and, make the N ^ oo 
limit while introducing the continuous variables x — i/N 
and y — j /N to arrive at 
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where ||a; — 2/|| represents the minimal distance on a circle 
of perimeter one. We can recognize the fractional integral 
/^^"and consequently write 

C(x) = i^r(l-a)/i-"(cosq(x)) . (9) 

In this large size limit, the a— HMF dynamics implies 
studying the evolution of the scalar fields q{x, t) and 
p{x, t) which are ruled by the fractional (non-local) partial 
differential equations 



dq 

'Ft 

dp ^1 

at ~ 2"' 

where /i = 1 — a. It has then been shown in (7) that 
stationary states are solutions of 



:p(x,t) 

-r(/i) (- sm{q)I'' (cos q) + cos(g)/^ (sin q)) 



d" cos q 

V cos q — — ; = . 

dx°' 



(10) 



where the operator D" stands for the fractional derivative, 
and that actually this property was shared with non 
stationary QSS's. All these results were obtained for the 
model in its long range version, meaning when a is smaller 
than one. 
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Fig. 1. Magnetization vs energy, for a — 1.5. A phase 
transition is displayed with Ec ~ 0.6. The transition 
point seem to be different from the long-range one 
which is Ec = 0.75, but the qualitative behavior of 
curve is the same. 

2. WHAT HAPPENS FOR 2 > a > 1 

We recall that systems are considered long range when the 
two body interacting potential V(r) decays at the most as 
1/r" with a < d, where d stands for the dimension of 
the embedding space. Having only one degree of freedom 
d = 1 for the a— HMF model. In these regards, considering 
situations where 1 < a < 2 is actually studying short 
range models. In fact TV is finite so there is no need of 
system size renormalization of the couling constant for 
a > 1. However something is peculiar about this lattice 
model. Indeed when considering the dynamics of a long 
range system, we would expect the force to be ruled by a 
1/r^ decay with f3 = a + 1, which is not the case for the 
lattice model for which the decay exponent is unchanged 
{l3 — a). Moreover given the particular importance of the 
microscopic dynamics and possible ergodicity breaking, 
one could naturally raise the question if the long range 
nature of a system is not ruled by the dynamics, which 
would then imply a system to be long ranged if /3 < d+1, 
which for the a— HMF model would imply a < 2. Most of 
the previous analysis of the model has been performed for 
a < 1, and can not simply be extended to 1 < a < 2. A 
first numerical analysis is therefore necessary. 

2.1 Phase transitions 

One peculiarity of one dimensional systems, is that there 
should not bare any phase transition if the interaction is 
short ranged. A first numerical study of the magnetization 
versus density of energy is performed for a = 1.5 in Fig. 1. 
The numerical integration of the microscopic dynamics 
is performed using a simplectic scheme, (optimal fifth 
order see(14)), a typical time step used is St = 0.05, and 
the initial conditions are Gaussian distributed. The fast 
Fourier transforms are done with the fftw libraries. As 
can be seen in Fig. 1 a phase transition is displayed at 
a transition point Cc ~ 0.6 which is then different from 
the universal value obtained for a < 1 which is Ec = 0.75. 
Preliminary results show actually that the critical point 
depends on the value of a and results seem to show that 
it approaches e = for a = 2. 

This existence of a phase transition beyond the "classical" 
long range threshold in one dimension had already been 
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Fig. 2. Top: Magnetization curve vs time for a = 1.0, N = 
2^8 and e — 1.2. During the QSS the magnetization 
is oscillating for a time ti, then it relaxes down to it's 
equilibrium value in a time T2. 

Bottom: — for different values of A^, a = 1 and e = 
0.6. The two times are approximately proportional, 
so the knowledge of ti gives a good approximation for 
T2- This proportionality is respected for all a values in 
this paper, even if the proportionality constant may 
vary. 

noticed for the Ising model by Dyson in the sixties, but this 
feature seems to favor the dynamical definition of what a 
long range system ought be. However an important feature 
to assert this new definition would be to find as well quasi- 
stationary states in this region of a's. 

2.2 QSS lifetime 

In long range interacting systems, generally, the limit N — > 
oo and i — > oo doesn't commute, so the thermodynamic 
limit is not unequivocally defined and one may end up in 
different equilibrium states depending on the order of the 
previous limits. Physically it's more feasible to compute 
the continuous limit before the time limit, so if the lifetime 
T of the QSS diverges with N it becomes the effective 
real equilibrium of the system, which in general is not 
obeying Boltzmann's statistics (4). We studied how the 
QSS lifetime scales with the exponent a in the decay 
parameter of the potential. 

First we studied the behavior of the lifetime r around the 
crucial value a = 1 to better understand the transition 
between a long range system and a supposedly short 
range one, but as mentioned in the case of a-HMF there 
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Fig. 3. Poincare sections for a == 1, e = 1.2 and TV = 2^°. 
Top refers to the QSS, while bottom represent the 
relaxation state between ti and T2. It can be easily 
seen that the second small isle disappear during the 
relaxation and the phase space becomes symmetric in 
q, thus ending the oscillations of the magnetization. 

are convincing arguments that the requirement for the 
emergence of long range behavior could be relaxed, and 
we can expect some long range feature to survive above 
a = 1. 

We considered the initial condition already used in (7) 
which was giving rise to a QSS, namely a long lived 
magnetized state above the critical energy. The initial 
condition used is all qi — 0, wile the p^'s are Gaussian. 
To characterize the QSS lifetime we monitor the behavior 
of the global macroscopic parameter magnetization which 
can be M — \J2i^^'^^\- ^^ Fig. 2 we show the behavior 
for a = 1, which is qualitatively representative for the all 
studied values of a in this paper. Here a first transition 
at t = Ti can be identified until which the system 
oscillates around an almost constant magnetization value, 
and beyond which the system starts to relax towards the 
equilibrium M = value. A second transition at i = T2 
is as well identified, it corresponds to the time at which 
the system finally reaches it's equilibrium state. As can be 
seen in Fig. 2, we find that these two values are almost 
linearly proportional for each a-value that we took into 
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Fig. 4. Scaling of ti with A^ for different a values and 
e = 1.2. Top shows the scaling for a — 0.9; middle 
refers to a — 1.0 and shows both ri and r2; while 
bottom refers to a = 1.1 . All curves seems to grow 
at least logarithmically. 

consideration, so we will refer to the first lifetime ri as the 
lifetime of the QSS, since is an order of magnitude faster 
to compute and we are interested only in the qualitative 
form of the scaling law for the lifetimes. 

The difference of the two dynamical regimes defined by 
the above thresholds are better understood when looking 
at the Poincare section captured in each of this regimes 
displayed in Fig. 3. At first the system is forms two 
distinct islands in the phase space, which start moving 
around and create the oscillations in the magnetization 
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Fig. 5. Magnetization curve vs time for a — 1.5, e — 1.0 
and increasing N values. The initial oscillation 
lifetime ri is independent of N and the system 
relaxes to equilibrium in a time which should be 
constant in the continuous limit. 



that characterize the QSS, and then one of the isles 
disappear during the relaxation period, while the phase 
space becomes symmetric in q thus ending the oscillations. 

Now we analyze the lifetime of the QSS versus the size 
of the system around the classical long range threshold 
a = 1, if the system is long range it should diverge with 
TV. The cases a = 0.9, a = 1 and a — 1.1 are displayed 
in Fig. 4. We can see that around a — 1 the scaling of t 
with TV approaches a logarithmic curve, meaning that the 
QSS survive at least until a — \ and maybe beyond as it 
appears as well true for a — 1.1. However when looking at 
the data for larger value of a, namely a = 1.5 the scaling 
of r appears to saturate as displayed in Fig. 5, where the 
magnetization curves obtained for a = 1.5 appear to all 
be identical no matter the size of the system. At this value 
of a there is still a short initial oscillation in M typical of 
the QSS, but now it's lifetime seems to be independent of 
N and finite so the system will reach the equilibrium state 
in a large enough time. Conversely we may expect that we 
may actually observe the same feature for a = 1.1, but 
that we have not seen yet the saturation in the lifetime as 
we were not able to simulate systems that would be large 
enough. In other words, we have as now not enough data to 
identify if there is a transition value between 1.1 < a < 1.5 
where the system becomes suddenly short-range dropping 
the logarithmic law of ri, or it still saturates at some 
larger value oi N > 2^° even for a — 1.1. Preliminary 
results shows that this saturation becomes quite fast for 
a ^ 1.2, 1.3 where even for low values of N it may still be 
possible to observe a scaling which is sub-logarithmic, but 
even if it did this would be only relevant for systems with 
an astronomical scale of constituants. 

However even if the QSS lifetime seems to saturate at some 
point, so that the system dynamics change into a short 
range one, the phase transition from a magnetized to an 
homogeneous state, typical of the long range regime, is still 
present (figure 1). Hence is appears that the dynamical 
definition may be more relevant for macroscopic features. 



such as the presence of a QSS (even with a finite Hfetime) 
or a phase transition, while the more classical statistical 
definition of a long range system corresponds actually 
to different dynamical behavior of the system and the 
existence of QSS with diverging lifetimes. 
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